Quantum-to-Classical Transition in Cavity Quantum Electrodynamics 
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The quantum properties of electromagnetic, mechanical or other harmonic oscillators can be re- 
vealed by investigating their strong coherent coupling to a single quantum two level system in an 
approach known as cavity quantum electrodynamics (QED). At temperatures much lower than the 
characteristic energy level spacing the observation of vacuum Rabi oscillations or mode splittings 
with one or a few quanta asserts the quantum nature of the oscillator. Here, we study how the clas- 
sical response of a cavity QED system emerges from the quantum one when its thermal occupation 
- or effective temperature - is raised gradually over 5 orders of magnitude. In this way we explore 
in detail the continuous quantum-to-classical crossover and demonstrate how to extract effective 
cavity field temperatures from both spectroscopic and time-resolved vacuum Rabi measurements. 

PACS numbers: 42.50.Pq, 03.67.Lx, 07.20. Dt, 44.40.+a 



Cavity QED [1 enables the study of the nature of 
matter light interactions in exquisite detail. It realizes 
an open quantum system in which the coupling to the 
environment is highly controllable. In a circuit real- 
ization of cavity QED [2], we carefully investigate the 
quantum-to-classical transition of a harmonic oscillator 
strongly coupled to a two level system by increasing the 
effective oscillator temperature. From measured vacuum 
Rabi splitting spectra and from time-resolved vacuum 
Rabi oscillations we consistently extract effective cavity 
field temperatures between 100 mK and a few K using a 
quantum master equation model as suggested in Ref. [3] ■ 
The dissipative quantum-to-classical crossover of a field 
mode coupled to a qubit was also studied theoretically in 
Ref. [4]. The emergence of classical physics from quan- 
tum mechanics and the role of decoherence in this process 
is an important subject of current research [5]. 

In our experiments a superconducting artificial atom 
is strongly coupled to a coplanar transmission line res- 
onator to realize a circuit QED system [2]. The trans- 
mon qubit [6; consists of two superconducting aluminum 
thin films weakly connected by two Josephson tunnel 
junctions. Its ground \g) and excited state |e) differ 
by the specific charge configuration of the two super- 
conducting islands. The transition frequency ^ gjC ~ 
^&EqEj{<&) — Ec is determined spectroscopically. Here, 
E(j/h sa 0.502 GHz is the single electron charging en- 
ergy, = i?j !max | cos (7T$/$ )I the flux controlled 
Josephson energy with Ej tinax /h ~ 14.4 GHz and $o is 
the superconducting flux quantum. The cavity has a bare 
resonance frequency of v v sa 6.44 GHz and a coupling lim- 
ited photon decay rate of k/(2it) 3.2 MHz. Optical 
microscope images of the sample are shown in Fig. [lja) . 

In our experimental setup [Fig. [ljb)] the coupled 
qubit/cavity system is prepared in its ground state \g,0) 
[Fig. [2ja)], with close to n — photons in the resonator 
by cooling the sample to below 20 mK in a dilution re- 
frigerator. Instead of increasing the physical tempera- 
ture of the sample to control the thermal occupation of 
the cavity, we apply thermal radiation only at its input. 



We approximate the one-dimensional Planck spectrum 
of the thermal field as constant with a power spectral 
density S n in the small relevant bandwidth around v r . 
S n is controlled by applying broadband microwave fre- 
quency white noise of controlled amplitude. In this case, 
we can assign an intra-cavity thermal photon number 
fhh — [exp (fiuj-r I 'kftT c ) — 1] _1 and thus an equivalent cav- 
ity field temperature T c to the externally applied white 
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FIG. 1: (a), Coplanar microwave resonator with two qubit 
slots and flux bias lines (top) and a single embedded qubit 
(violet) shown on an enlarged scale (bottom), (b), Circuit 
diagram of setup. Input (left) is at room temperature (RT) 
with microwave sources for measurement (i/ plo b c ), qubit spec- 
troscopy and qubit drive (^ spC c), variable attenuation (D) 
quasi thermal field source (Swhitc) at v Y and an arbitrary 
waveform generator for fast flux biasing (Idc) through the 
on-chip flux line (B). At 20 mK the qubit is coupled via a 
capacitance C g to the transmission line resonator between the 
capacitances Ci n and Cout- The transmitted microwave tone 
is amplified, down-converted with a local oscillator (LO) and 
digitized (ADC). 
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noise. Owing to the high internal quality factor of the res- 
onator, the field does not thermalize on the chip, which 
allows us to control the effective temperature of the res- 
onator field to up to T c ~ fOOK. In addition our setup 
allows for the phase sensitive detection of the quadrature 
amplitudes of a weak coherent probe tone that popu- 
lates the resonator with n pro b ^0.1 photons on average 
while effectively rejecting all uncorrelated thermal radi- 
ation from the detection system. 

At the lowest measured cavity field temperature 
T c ~ 100 mK, a clear vacuum Rabi mode splitting 
is observed in a linear response cavity transmission 
measurement, see Fig. [2^c) (dark blue lines). Two 
Lorentzian lines characteristic for the dressed states 
| n, ±) = (\g, n) ± |e, n — 1}) /v2 with excitation number 
n = 1 separated by twice the dipole coupling strength 
g gj e/(2ir) = 54 MHz are clearly observed at the fre- 
quencies t'go.ii- This indicates that the strong cou- 
pling regime, g gie ^> K, 7, with the qubit relaxation rate 
7/(2tt) « 0.6 MHz, is realized. 

While the observation of the simple vacuum Rabi mode 
splitting for n — 1 could in principle be interpreted as a 
normal mode splitting using a semiclassical model, the 
observation of additional transitions to higher excited 
dressed states with coupling increased by a factor \fn 0- 
[TD] is an unambiguous signature for the quantum nature 
of the system [LTHT5] . In our experiment these transi- 
tions are observed at frequencies v\±p±, see the level di- 
agram in Fig. ffia ) and measurements in Fig. [2]^c) (green 
lines) [and Fig.j3]ja) (blue and green lines)], when slightly 
increasing the effective cavity temperature and thus pop- 
ulating excited states of the Jaynes-Cummings ladder. In 
this case the weak coherent probe tone excites transitions 
between higher levels. 

At high temperatures the Jaynes-Cummings system is 
excited to high quantum numbers and many transitions 
can be accessed by the weak probe tone [3J [TTHT4"] . The 
resulting linear response transmission spectrum can be 
understood as a sum of individual level to level transi- 
tions that overlap. At a mean thermal photon number of 
n t h > 1 the transitions from the ground state f g o,i± are 
almost saturated and therefore only weakly contribute 
to the observed cavity transmission spectrum (Fig. [2^c) , 
yellow lines). At even larger effective temperatures all 
available transition frequencies are densely spaced close 
to the resonator frequency, see level diagram in Fig. J2jb) 
and measurement in Fig. [2j(c) (red lines) . At these tem- 
peratures all transitions which would show the nonlinear 
\Jn signature of the single qubit-photon interaction are 
saturated and the transmission spectrum does not carry 
any information about the intrinsic quantum nature of 
the system. 

The resulting harmonic spectrum resembles that of two 
uncoupled linear oscillators, a situation where not the 
absolute energies but the energy differences of the sys- 
tem are degenerate [3j [15] . In contrast to the regime of 
strong coherent driving, where nonclassicality can persist 
at large excitation numbers [TCJ [T7] , the phase noise of 
the thermal field used in our experiments renders the two 
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FIG. 2: (a), Energy level diagram of dipole coupled dressed 
states \n, ±) for excitation numbers n = 0, 1 and 2 (black 
lines), uncoupled qubit and cavity energies (grey dashed lines) 
and allowed transitions at frequencies v s o,i± and v\± } i± (blue 
arrows), (b), Dressed state diagram for large excitation num- 
bers n > 280 and allowed transitions at the resonator fre- 
quency u r (red arrows), (c), Measured cavity transmission 
yl 2 / J 4n lax for intracavity thermal photon numbers 0.05 (blue) 



< 



nth 



< 



400 (red) and a fit to a Lorentzian line (black). 
Measured data sets are normalized to the transmission am- 
plitude j4 max obtained when the qubit is maximally detuned. 
The data is offset and colored for better visibility, see inset in 
Fig. [4] for color code. 



high-excitation dressed-state ladders indistinguishable. 
As a result at high temperatures the spectrum is that of a 
classical resonator and we find excellent agreement with 
a Lorentzian line fit of width ft e g/(27r) = 4 MHz close to 
the intrinsic resonator line width (Fig.^c), black line). 

A quantitative understanding of the measured results 
is obtained by numerically solving the Markovian master 
equation 



C[p] =--[U,p] + (nth + l)nV[&]p 
n 
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where the thermal photon number n t h or equivalently the 
cavity field temperature T c is extracted as the only fit pa- 
rameter. Here the coherent dynamics is described by the 
Jaynes-Cummings Hamiltonian % with I transmon lev- 
els, see [B], and without a drive term, which is justified 
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FIG. 3: (a) Measured (dots) and calculated (lines) cavity transmission (same data as in Fig. 2(c)) shown for applied S n = 
— 221.5 dBm/Hz (bottom) to —190 dBm/Hz (top) in steps of 4.5 dBm/Hz. Extracted thermal photon numbers nth and cavity 
temperatures T c are indicated, (b) Measured qubit excited state population P c as a function of the resonant cavity interaction 
time t (dots) and master equation simulation (lines) for S n = —214 dBm/Hz to —202 dBm/Hz in steps of 3 dBm/Hz. (c), 
Similar measurement as in (b) with qubit prepared in the excited state for S n = —214 dBm/Hz. Inset: Vacuum Rabi pulse 
sequence. 



in the linear response limit. The three damping terms 
model the loss of cavity photons at rate k, the intrin- 
sic relaxation of the transmon excited state |e) at rate 7 
and the creation of cavity photons due to a thermal bath 
with occupation n t h- 2? is the usual Lindblad damping 
superoperator. Note that neither qubit dephasing nor an 
independent qubit thermal bath has significant influence 
on the agreement between the theory and the presented 
spectroscopic results. Therefore, these two terms have 
been omitted. Equation ([I]) is solved by exact diago- 
nalization and the transmission spectrum is calculated 
similar to Ref. [3]. 

We observe that the measured and calculated spec- 
tra agree very well over a large range of applied thermal 
noise powers (Fig. [3]). This allows to extract the effec- 
tive cavity field temperature (Fig. |4|. For mean thermal 
photon numbers n t h ^ 30 the accuracy of our numeri- 
cal calculations is limited by the finite size of the Hilbert 
space spanned by 200 resonator states and 6 transmon 
states. In the large photon number limit the qubit only 
negligibly perturbs the empty cavity Lorentzian spec- 
trum [Fig.[2[c)], because the cavity and cavity dissipation 
terms in Eq. ([!]) scale linearly with n, while the qubit 
photon coupling term scales only with yjn, see also |18j . 
In the limit where the coupling is smaller than the damp- 
ing y/ng g ^ c < tlk the dominant terms in Eq. ([!]) describe 
a damped harmonic oscillator with a Lorentzian spec- 
trum. A classical harmonic oscillator is furthermore a 
good description of our cavity QED system if the rel- 



ative y/n nonlinearity, i.e. the frequency shift due to 
the qubit, is smaller than the relevant dissipation rate 
2 <7 gj e(V n + 1 — v 7 ") < K - F° r large photon numbers 
these two criteria arc equivalent and in our experiment, 
the classical limit is entered for mean thermal photon 
numbers n t h > (<? g , c / K ) 2 ~ 280. 

We have also performed time domain vacuum Rabi os- 
cillation measurements in the presence of thermal pho- 
tons. In these experiments, we first detune the qubit from 
the resonator by 0.5 GHz where the applied quasi ther- 
mal noise is filtered and does not affect the qubit state. 
We then apply an amplitude shaped current pulse i^c via 
an on-chip flux line, see sample and setup in Fig. [TJ to 
tune the qubit in resonance with the cavity field for a 
variable time r, see pulse sequence in inset of Fig. [3^c) . 
The qubit excited state population P c is then determined 
in a dispersive measurement [19] . 

For the lowest cavity temperatures we observe low con- 
trast, long coherence vacuum Rabi oscillations induced 
by the weak thermal field only [Fig. [3^b)] . When the 
number of thermal photons is increased the amplitude of 
the coherent oscillations increases while their coherence 
time decreases, which is expected for a thermal distribu- 
tion of photon numbers. At long interaction times r, the 
qubit population shows a temperature dependent satu- 
ration that approaches the value of 0.5 when the qubit 
is strongly driven by the thermal field at high tempera- 
tures. We have repeated the experiment preparing the 
qubit initially in the excited state using a 7r-pulse and 
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FIG. 4: Temperature T c and photon number nth extracted 
from measured Rabi spectra (full dots) and Rabi oscillations 
(crosses) versus noise power spectral density S n or equiva- 
lently attenuation D. Theory (black line). Inset: extrapola- 
tion to higher T c provides temperature scale (color coded) for 
the measurements shown in Figs.[2^c) and [3] 

measuring vacuum Rabi oscillations at low applied noise 
power [Fig. [3jc)] . In this case high contrast oscillations 
are observed albeit with a small offset in the qubit popu- 
lation visible at large r caused by the thermal field. From 
these data the cavity temperature has been extracted by 
solving a time dependent master equation, see solid lines 
in Figs. (3^b) and (c). In this case 2 qubit states and 
up to 6 cavity states are considered and qubit dephasing 
is included. The deviations at short times r are due to 
imperfections in the qubit tuning pulses. 

The cavity field temperatures T c extracted from both 
spectroscopic (full dots) and time-resolved measurements 
(crosses) plotted on a logarithmic scale versus the applied 
noise spectral density S n are in close agreement (Fig.Ej). 
Both data sets are explained consistently considering 
that the cavity thermal photon number n t h — S n /fkJ+rio 
is the sum of the applied quasi thermal noise S n /Kui ex- 
pressed as a photon number and a small thermal back- 
ground field of n = 0.04 photons (T = 95 mK). The 
measured intracavity photon number n t h is consistent 
with the noise applied to the input of our setup and the 
expected attenuation of the line and also with calibration 
measurements based on the known qubit ac-Stark shift 
|20j . The background photon number is set by incom- 
plete thermalization of the input line and can easily be 
decreased by adding additional cold attenuation. 

The inset of Fig. [4] shows an extrapolation of the tem- 
perature to the largest applied noise powers and provides 
a calibration for the measurements shown in Fig. J2jc) 
where the same color code is used (also used in Figs. [3] 
and [4| . Accordingly, the measurements at the highest 
temperatures were conducted with approximately ~ 370 
thermal photons in the cavity. This corresponds to an 
effective cavity temperature of T c ~ 115 K. We clearly 
observe the transition to a classical harmonic resonator 
at n t h ~ 370 > (gg^/w) 2 « 280 as indicated by the 
Lorentzian spectrum [Figj2|c)]. 



We have demonstrated a quantitative understanding of 
the transition from the quantum to the classical response 
of a cavity QED system. Moreover, we have presented a 
viable approach to measure the temperature of an elec- 
tromagnetic field in a cavity over a wide range. A differ- 
ent approach to extract the resonator temperature in the 
dispersive regime of circuit QED has been proposed [5T] 
and measurements of the qubit temperature have been 
demonstrated [55]. Related experiments have been per- 
formed using a semiconductor cavity QED system [23] . 
In these systems, entanglement and decoherence at ele- 
vated temperatures can be studied in future experiments 
in the context of quantum information. 
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